Quasiparticle Transport in the Vortex State of Unconventional Superconductors 
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I consider the problem of the vortex contribution to quasiparticle transport in unconventional super- 
conductors with line nodes, and argue that the magnetic field dependence of transport coefficients 
is fixed by the same scattering processes which limit low-temperature transport. I give no detailed 
calculations, but show that qualitatively correct results may be obtained in the limit of low tem- 
peratures and fields by simple physical arguments, based on estimates of the density of states and 
relaxation time in analogy with the zero-field case. I conclude with a brief discussion of the problem 
of anisotropy of the field dependence and influence of experimental geometry. 



PACS Numbers: 74.25.Fy, 74.72.-h,74.25.Jb 

Introduction. For some time, transport measurements 
in the Abrikosov vortex state of cuprate superconduc- 
tors have exhibited magnetic field and temperature de- 
pendences which seemed difficult to incorporate into the 
standard picture developed for classic superconductors. 
|E| In this model, if vortices are pinned, quasiparti- 
cles carry all electrical current and entropy, and their 
transport is limited exclusively by scattering from vortex 
cores. The cores play a crucial role because of a rela- 
tively high density of bound states created there by the 
local order parameter supression, and because extended 
quasiparticle states are completely gapped. 

In the past three or four years, a consensus has been 
established that the order parameter in the hole-doped 
cuprate superconductors has <i-wave symmetry. [0] It is 
therefore natural to ask if the new symmetry class, and 
in particular the existence of gap nodes on the Fermi sur- 
face, can invalidate the "vortex scattering" model used 
to analyze quasiparticle transport in classic superconduc- 
tors. In this paper, a summary of two lectures given at 
the February 1998 Workshop on Strongly Correlated Sys- 
tems of the Asia Pacific Center for Theoretical Physics, 
I argue that this is in fact the case: that the existence of 
low-energy quasiparticle states due to the order parame- 
ter nodes allows one to neglect vortex scattering at small 
temperatures and fields, and focus instead on the rear- 
rangement of nodal quasiparticle occupation numbers oc- 
casioned by the Doppler shifts of quasiparticle energies 
in the superflow field of the vortex lattice. Relevant scat- 
tering processes are the same which dominate zero-field 
transport, leading to the possibility that the successful 
"dirty d-wave" theory of the H = superconducting 
state transport in HTSC can be extended to the vor- 
tex state without the introduction of further parameters. 
Predictions of the theory for the field dependence of, e.g., 
the thermal conductivity in the vortex state, are consis- 
tent with measurements made to date at higher temper- 
atures. At low temperatures, where the predictions are 
drastically different from those of the vortex scattering 
model, preliminary experimental results support the the- 
ory, but more experiments are needed. 

I begin by providing the qualitative background neces- 
sary to understand the temperature dependence of H = 



transport coefficients in the cuprates. I then summarize 
the arguments why extended quasiparticle states may be 
assumed to dominate the vortex bound states, and dis- 
cuss the semiclassical approximations used to treat them. 
I show that at T = 0, transport coefficients increase 
quasilinearly with field, but at finite temperatures ex- 
hibit a nonmonotonic field dependence. For the most 
part I restrict myself in this paper to qualitative argu- 
ments which can be shown to provide the same results as 
obtained by a more complete treatment; in doing so I am 
mostly summarizing work with C. Kubert which has ap- 
peared elsewhere. ||-||] However, in section 4, I consider 
the influence of relative configuration of magnetic field 
and current, and give as a concrete example new results 
for thermal transport as a function of field in the basal 
plane of a d x 2_ v 2 superconductor. 



I. INTRODUCTION: ZERO-FIELD TRANSPORT 

At sufficiently low temperatures, the quasiparticle 
transport time will be determined by impurity scattering. 
If one restricts one's attention to the s-wave scattering 
channel, but assumes an unconventional superconducting 
state in which = 0, Pethick and Pines || showed 

in the context of heavy fermion superconductivity that 
1/t(w) varies as N(u>) in the conventional weak scatter- 
ing (Born) limit, and as iV(w) -1 in the strong scattering 
(unitarity) limit. The strength of individual scattering 
events therefore has profound consequences for transport 
coefficients L, which may be estimated as 



N{oj)t(uj) ~ N(T)t{T) 



(1) 



N(T) 2 unitary 
Ln — const. Born 



where Ln — L(T C ). 

Experimentally, in both heavy fermion and cuprate su- 
perconductors, it is found that in good samples, bulk 
transport coefficients are orders of magnitude smaller 
than their normal state counterparts, thus leading to the 
(hitherto microscopically unjustified) ansatz of unitarity 
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scattering by simple defects in these sy terns. Formally, 
such calculations are carried out either in the kinetic 
equation approach of Pethick and co-workers, R] or in the 
self-consistent i-matrix approximation (SCTM). ||[| In 
the unitarity limit, both approaches yield similar results 
for temperatures and energies above a scale 7 ~ \J AoT, 
where Ao is the gap maximum and T is the impurity 
scattering rate in the normal state. At low tempera- 
tures, only the SCTM yields true gapless behavior (finite 
residual density of states 7V(0)), as well as the remarkable 
prediction of universal transport, |!o| ], recently confirmed 
by experiment in thermal conductivity measurements on 
Zn-doped YBCO. Briefly, this occurs because zero- 
energy quasiparticle states acquire a residual broaden- 
ing 1/t(0) = 7 in the presence of impurities due to the 
nodes of the gap, and the residual density of states is 
also found to scale with 7/A0. Thus L is generically of 
order No(-f / Ao) / '7 ~ Nq/Aq, i.e. L is independent of 
the impurity concentration to leading order in the dis- 
order r. To give an explicit example, in a d x 2_ y 2 state, 
Afc = Ao cos 2cf> over a model cylindrical Fermi surface, 
the thermal conductivity in zero applied magnetic field 
for heat current jg in the plane is approximately 



K el {T) 



7nT 2 

5A r 



In 



1.14A 



7T J 

T 



(2) 



where kqq/T = ttNqVp/(6Aq) is the universal thermal 
conductivity. [£^]lj] 

At higher temperatures in the cuprates, the mean free 
path due to inelastic scattering becomes smaller than 
that due to impurity scattering. DC transport coeffi- 
cients generically then show the following structure: a 
rise at low temperatures from the universal value to a 
maximum at intermediate temperatures, followed by a 
decrease to the normal state value at T c , typically an 
order of magnitude or so less than the peak in clean sys- 
tems, but an order of magnitude greater than the coef- 
ficients at the lowest temperatures. (See Fig. 1) The 
existence of the maximum has been qualitatively under- 
stood for some time as arising from the competition be- 
tween a collapsing inelastic scattering rate l/ri ne i be- 
low T c due to an opening of a gap in the spectrum of 
excitations responsible for inelastic scattering, and the 
decrease of the number of quasiparticle carriers at low 
temperatures. |l5| , |l6| ] At low temperatures there may be 
problems with the details of this simple argument, due 
to the energy dependence of the scattering rate due to 
impurities, l/Ti mp (as discussed above), but the overall 
features are well-described by models which add impu- 
rity and inelastic scattering rates incoherently and in- 
clude the collapse of the inelastic scattering below T c in 
some form. My co-workers and I argued in [[l7| that at 
sufficiently low temperatures, when quasiparticle states 
are well-defined, l/ri„ e ; should vary as T 3 for the simple 
reason that electron-electron collisions in a normal Fermi 
liquid with constant particle number leads to T 2 , and 
there is one additional factor of the temperature arising 



from the linear d-wave density of states. This estimate is 
borne out in a weak-coupling Hubbard model calculation 
of 1/t, jl8| which yielded excellent results in compari- 
son to both low-frequency microwave conductivity |l7| l 
and thermal conductivity [ fl9[ measurements, including 
disorder dependence, over the entire temperature range. 
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FIG. 1. Schematic behavior of DC transport coefficients in 
the cuprates. 



II. QUASIPARTICLES IN THE 
UNCONVENTIONAL VORTEX STATE 

I now turn to the question of transport in the vortex 
state, focussing for consistency on longitudinal thermal 
conductivity measurements, although microwave conduc- 
tivity measurements should share many of the same fea- 
tures. Thermal conductivity as a probe of the quasipar- 
ticle system has the advantage that in most cases of cur- 
rent interest vortex motion may be ignored. While the 
phonon contribution is often very difficult to determine, 
it seems likely that it is field independent, so the quasi- 
particle system is directly probed by measurement of the 
field dependent part 6n(T; H). This question is still open, 
however. There are several measurements of n at temper- 
atures down to 10K or 20K in YBCO-123 and BSCCO- 
2212, a typical recent one being that of Freimuth, |2C| ] 
which shows at fixed temperture T ~ 0.1T C that the 
thermal conductivity decreases with application of a field 
with sublinear power or log behavior. The qualitative 
behavior at higher temperatures is not inconsistent with 
the conventional vortex scattering model, in the sense 
that the quasiparticle mean free path shortens with in- 
creasing field as might be expected if the vortices are 
moving closer together. However, the form of the field 
dependence is unusual, since, within the picture just al- 
luded to, the number of scatterers should scale linearly 
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with H. Furthermore, recent measurements at lower T 
on both YBCO (2l] and UPt 3 , |f] a heavy fermion sys- 
tem also with an order parameter with line nodes, have 
exhibited a thermal conductivity increasing quasilinearly 
with field. This can never be explained within a vortex 
scattering picture. 

Vortex scattering in the cuprates is in fact likely to 
be negligible at low fields for several reasons. The den- 
sity of bound states in a state with gap nodes was found 
to be significantly smaller than the density of extended 
states by Volovik |^3| when averaged over the entire vor- 
tex. This calculation accounted for the gap nodes, but 
treated the bound states as a continuum and thus must 
be regarded as an upper bound to the true contribution. 
In fact, in the cuprates the bound state separation or 
"minigap" S ~ Aokp 1 /£o is very large compared to clas- 
sic superconductors, and it is unlikely at low T that more 
than a single state/vortex is occupied. The net cross sec- 
tion for quasiparticle scattering is therefore likely to be 
very much smaller than in the classic case. Finally, I note 
that at low temperatures the picture of independent scat- 
tering events by vortices must break down in the perfect 
vortex lattice. The quasiparticle-vortex mean free path 
^vortex will then arise solely from the disorder-induced or 
residual thermally excited fluctuations of the ideal peri- 
odic Abrikosov lattice. At H — IT, the intervortex dis- 
tance of 2R = ^(^y^a'^H^/H) 1 / 2 , which provides 
an absolute lower bound for £ vor tex, is about 25 times the 
coherence length £o, or perhaps 400 A in HTSC. While 
in YBCO the mean free path due to defect scattering 
and electron-electron collisions ltr at low T may be as 
long as 3000A, in BSCCO ltr is a factor of ten or so 
smaller, p4| and therefore probably smaller than l vor tex 
over most of the field range of 10T or so in current experi- 
ments. Note that this same argument applies to Andreev 
scattering by the vortex flow field as considered in Ref. 
pB| , since again vortex disorder, rather than the inter- 
vortex distance, sets the mean free path. 

// mean free paths due to vortex scattering are indeed 
much longer than in classic superconductors, other scat- 
tering processes, such as impurities and electron-electron 
collisions, will limit the mean free path at low fields 
H <C H C 2- I will simply assume this for the discussion be- 
low, and beg the question of the range of fields for which 
such an assumption may be justified. 

In his seminal paper, Volovik showed that the domi- 
nant contribution of the extended quasiparticle states to 
the density of states took the form 



N(0;H)~N J — 



H 

H C 2 



(3) 



Such a term was indeed extracted from specific heat data 
on YBCO by Moler et al., j26) at the time a crucial piece 
of evidence in support of the identification of d-wave sym- 
metry in YBCO. This result can be reproduced easily if it 
is recalled that 1) the density of states in the pure d-wave 



state varies as N(u>) ~ w/A at low energies lj <C A ; 
and 2) that a quasiparticle of momentum k moving in 
the superflow field of a single vortex, v s (r) ~ h/(2mr)a 
(a is the azimuthal angle in real space, e.g. in the plane 
if H || c) experiences a Doppler shift u> — > u> — v s • k in 
the lab frame. We then have at uj = 



N(0;H) 



A 



H 



1 

R 2 



dr i 



(4) 



where in the second approximate equality an average 
(...) h has been performed over the vortex unit cell, of size 
R. Since the integrand is constant, N(Q;H) ~ R^ 1 ~ 
\f~H follows immediately. 



III. THERMAL CONDUCTIVITY AT H > 

Similar arguments can be applied || to the kernel of 
the thermal conductivity Kubo formula. At T = 0, anal- 
ogous to the discussion in Sec. 1, we have 



N(0;v s )~v s ; r(v s ) ~ N(0; v.) ■ 



(5) 



since the integrand now varies as 1/r rather than as a 
constant. Note in this case the average requires a lower 
cutoff, of order £o, which reflects our ignorance of how 
the extended quasiparticle wavefunctions vanish in the 
core region. Furthermore, we implicitly assumed in Eq. 
(5) that the local conductivity n(r) depended only on r: 
this is true only for the special case jq || H. A treatment 
of the opposite case, jg 1 H shows that the measured 
(averaged) k is much smaller since once necessarily adds 
thermal resistances in series rather than in parallel in this 
case. H A full numerical evaluation shows a behavior 
slightly more linear in field than given by Eq. (5), due to 
higher contributions of order v*, which give rise to true 
linear terms when averaged. 

At higher temperatures, terms decreasing in field arise 
due to 1) logarithmic energy dependence of the unitar- 
ity limit scattering rate (neglected in Eq. (1)), and 2) 
inelastic scattering. The competition between increasing 
and decreasing contributions leads in this case to a min- 
imum in the field dependence, at H ~ T 2 in the case 
of 1). Determination of the exact form including inelas- 
tic scattering has not yet been worked out in detail, but 
we know these effects are required for the understanding 
of cuprates in zero field above about IK or so. Simple 
estimates show that the minimum in this more realistic 
case is shifted out to fields higher than heretofore mea- 
sured; one can thus via this method understand both the 
magnitude and form of the (apparently) monotonically 
decreasing field dependences reported in, e.g. Ref. p0|. 
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FIG. 2. Schematic behavior of k(T; H)/T at various tem- 
peratures for a system with order parameter line nodes. Solid 
lines: clean system at T — 0, < T < T c ; dashed line: dirty 
system at T = 0. 



IV. ANISOTROPY 

In the above estimates, I have supressed consideration 
of spatial angular degrees of freedom in order to extract 
the qualitative field dependence. In general, the four- 
fold structure of the gap in momentum space fixed to 
the crystal field axes will induce an angular anisotropy 
in the current response due to the coupling v s • k. On 
the other hand, I have already mentioned in Sec. Ill a 
further source of anisotropy due to the different types of 
spatial averaging necessary for different directions of jq 
relative to H. By itself, this leads generically to a max- 
imum in the thermal conductivity when jq || H, consis- 
tent with observations by Yu et al. |25| and Aubin et al. 
f27f . Unfortunately, this qualitative feature is also com- 
mon to vortex scattering theories. It would be useful to 
consider an effect which is clearly due entirely to the gap 
anisotropy. 

Such an example is provided by taking the thermal 
current ]q parallel to the field H in the plane, with both 
making an angle e with respect to the rc-axis. In a co- 
ordinate system in which the quasiparticle momentum 
is k = (cos 0, sin 0, 0), the superfluid velocity is then 
v s = (— sine sin a, cos e sin a, cos a) (reminder: a is the 
vortex winding angle). The quasiparticle energy shift is 
then 



v s ■ k = v s k,F sin a sin(0 — e) 
v s k F 



(6) 



V2 



sina(± cos e ± sin e), 



where the final approximate equality is justified becasue 
the kernel in the Kubo formula is strongly peaked at the 
nodal momenta k n with <f> = 7r/4, 37r/4, .... Repeating 
the previous calculation, including the 
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FIG. 3. Anisotropy of [«(T; H) - k(T; 0)]/koo as function 
of field direction e (see insert) relative to the crystal axes for 
T = and T = 0.2T C and H/H c2 =0.01. Weak coupling value 
Ao/T c — 2.14 was assumed for the d x 2_ y 2 order parameter. 

angular average, and remembering to symmetrize over all 
nodes, one finds (up to factors of order unity) 



*(g) 
T 



k 



2 -271- 



(7) 



da sin 2 a max(cos 2 e, sin 2 e) 



n 



so the thermal conductivity is seen at T — to exhibit 
4-fold oscillations in the field dependence as a function 
of field angle in the plane with a mimimum in the nodal 
directions. In Figure 3, I plot a full evaluation of the low- 
energy approximation given in Ref. j^] as a function of 
angle in the basal plane. Note that a traditional vortex 
scattering model would predict no effect in this case. At 
a higher temperature, the model of Ref. Q predicts a 
decreasing n(H) for the field specified, and consequently 
maxima along the nodal directions; thus the oscillations 
appear similar to those observed in an experiment with 
slightly different geometry (but with H in the plane) by 
Yu et al. ]2q[ , and attributed to Andreev scattering. A 
comparison of these two effects is in progress. 



V. CONCLUSIONS 

I have tried to give a somewhat pedagogical summary 
of the ideas currently used to interpret transport exper- 
iments at low temperatures in d-wave and other uncon- 
ventional superconducting states. Back-of-the-envelope 
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estimates for the field dependence of thermal conductiv- 
ity in a field were shown to reproduce the more detailed 
calculations of recent work suggesting that the field de- 
pendence of transport coefficients in heavy fermion and 
cuprate superconductors can be well described by ignor- 
ing vortex scattering and including only H = scatter- 
ing processes, but shifting quasiparticle energies in the 
vortex superflow field. Fourfold oscillations of transport 
coefficients as a function of the angle between magnetic 
field direction H || j and the crystal axes in the plane 
of a d x 2_ y 2 superconductor were shown to be a natural 
consequence of this picture. 
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